Production of photons in a bouncing universe 



J. M. Salim 1 , S. E. Perez Bergliaffa 1 ' 2 and N. Souza 2 
Centro Brasileiro de Pesquisas Fisicas, Rua Xavier Sigaud, 
150, CEP 22290-180, Rio de Janeiro, Brazil, and 
2 Departamento de Fisica Teorica, Instituto de Fisica, 
Universidade do Estado de Rio de Janeiro, CEP 20550-013, Rio de Janeiro, Brazil 



Using a new non-singular solution, it is shown that the production of photons in dilaton electro- 
dynamics in a cosmological setting is increased if the effect of matter creation on the geometry is 
taken into account. This increment may be related to the problem of the origin of magnetic fields 
in the universe. 
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Introduction 



The origin, evolution, and structure of the magnetic 
fields present in galaxies and clusters of galaxies are 
amongst the most important open issues in astrophysics 
and cosmology. Typically, the magnetic field present in 
galaxies is of the order of a few /iG, and is coherent on the 
galactic scale. The observed field in the case of clusters is 
correlated over 10-100 kiloparsecs, and can be as high as 
tens of ^G. The standard mechanism to account for these 
fields is the dynamo 1], which amplifies small seed fields 
to the abovementioned values Q. The seeds that fuel 
the dynamo may have an astrophysical origin, or may be 
primordial. An astrophysical mechanism for the genera- 
tion of these pre-galactic magnetic fields is the Biermann 
battery [3j , which works for the generation of seeds to be 
amplified by the dynamo in galaxies, but can hardly ac- 
count for the fields in clusters 0, . Other astrophysical 
mechanisms (involving for instance starbursts, jet-lobe 
radio sources, or accretion disks of black holes in AGNs 
0) seem to need pre-existing magnetic fields. The cur- 
rent prevalent view is that the magnetic fields observed 
in galaxies and clusters have a primordial origin (see Q 
for a review). The processes that may account for this 
origin can be divided in two types: causal (those in which 
the seeds are produced at a given time inside the horizon, 
like QCD and EW phase transitions Q), and inflation- 
ary (where correlations are produced outside the horizon 
8J). In the latter, vacuum fluctuations of the electro- 
magnetic field are "stretched" by the evolution of the 
background geometry to super-horizon scales, and they 
could appear today as large-scale magnetic fields. How- 
ever, since Maxwell's equations are conformally invariant 
in the FRW background, the amplification of the vacuum 
fluctuations (which amounts to particle production) via 
inflation can work only if conformal invariance is broken 
at some stage of the evolution of the universe. There 
are several ways in which this invariance can be broken 
0: nom-minimal coupling between gravitation and the 
electromagnetic field 0, Hfl UM > quantum anomaly of the 
trace of the stress-energy tensor of electrodynamics 03 , 
coupling of the EM field to a charged scalar field H ^1 1 
exponential coupling between a scalar (whose potential 



drives inflation) and the EM field 01 ■> an d a non-zero 
mass for the photon Il5l. Yet another possibility is dila- 
ton electrodynamics |16( , in which there is a scalar field 
(the dilaton) which couples exponentially to an Abelian 
gauge field. In this model the inflationary expansion is 
driven and not by the potential but by the kinetic term 
of the scalar field. The exponential coupling is naturally 
implemented in the low-energy limit in string theory |17| , 
and in Weyl integrable spacetime (WIST) [lSf. and can 
be viewed also as a time dependence of the coupling con- 
stant, an idea considered first by Dirac |l9|. This avenue 
has been pursued by Giovannini in a series of articles 
[20|. More recently, a model with both a dilaton and 
an inflaton was analyzed in |21|. and later generalized to 
a noncommutative spacetime |23. In all these articles, 
different aspects of photon production have been ana- 
lyzed. We will focus in this contribution on two definite 
items. First , we shall address the creation of photons 
in a non-singular universe described by an exact solution 
of the field equations. More specifically, in this solution 
the passage through the bounce is described by the equa- 
tions of the model without resorting to unknown (Planck- 
scale) physics. The results of this calculation are to be 
compared with those of the second item, coming from 
a more complete model. We shall take into account in 
this second part the influence of the creation of matter 
in the squeezing of the vacuum state, a problem that 
seems to have received little attention in the literature. 
It will be shown that as a result of matter creation, there 
is an increment in the number of photons originating in 
the stretching of the vacuum fluctuations. Thus this ef- 
fect may be of importance in the generation of primor- 
dial magnetic fields. The analysis will be done in the 
framework of a phenomenological model, in which the 
creation of ultrarelativistic matter is restricted to an in- 
terval centered in the bounce. Two interesting features 
of the solution representing this model are that it goes 
automatically into a radiation regime after a short time, 
and it displays a constant value for the dilaton after en- 
tering the radiation era. We shall start in Sect.l by de- 
ducing a Hamiltonian for the EM field in a curved back- 
ground coupled to the scalar field. The quantization of 
this Hamiltonian system will be presented in Sect. 2, fol- 
lowing the method of the squeezed states. In Sect. 3 we 
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shall study the production of photons in the absence of 
matter (i.e. with only the scalar field present). Variation 
of the photon number in a model that takes into account 
matter creation will be discussed in Sect. 4. We close with 
a discussion in Sect. 5. 



In the comoving frame, the EOM along with Eqn.QJ be- 
come 



D % = 0, 



(7) 



I. FIELD EQUATIONS 



B a \\ a = 0, 



(8) 



A time-dependent gauge coupling is a generic feature 
of 4-dimensional theories obtained by compactification of 
some theory formulated in a spacetime with more than 4 
dimensions, such as Kaluza-Klein theory [2^ and string 
theory j24[. Time-dependent gauge couplings are also 
present in WIST [3. In all these cases, the action can 
be conveniently written in the form 



(1) 



where u> is the the dilaton (in the case of string theory) 
or the scalar field associated to Weyl geometry (in the 
case of WIST) and is an Abelian field. The function 
f(w) will be left unspecified for the time being, but later 
we will set f(uj) = c~ 2uJ , which corresponds to the case of 
string theory and WIST. The equations of motion (EOM) 
that follow from this action are: 



(f(cj)Fn\\ v =o, 



F a ?, = 

11/9 ' 



(2) 



(3) 



where the derivative wrt the background geometry is de- 
noted by " 1 1" , and the dual of the electromagnetic tensor 



F is defined as 



F af} = -r\ a ^ v F^. 



(h«D a ) + -9D x 



h?B a 



\0B X 



0. 



(9) 



(10) 



with 9 = 3d/ a, a the scale factor of the FRW metric, a 
dot represents derivative wrt cosmological time, and h" 
is the metric of the 3-space orthogonal to v^. As shown 
in the pioneering work by Lifshitz pl| , it is useful to ex- 
pand the solutions of linear differential equations defined 
on a curved manifold in the corresponding spherical har- 
monics basis. In order to do so we will first split the 
equations in space and time parts. The vectors describ- 
ing the electromagnetic field are already space vectors, 
since E a V a = D a V a = B a V a = H a V a = 0. From now 
on we shall work in conformal time, with the FRW metric 
given by 

ds 2 = a(ri) 2 [dr] 2 — 7y (x)dx l dx 3 ] , 

where 7^ is the metric of the hypersurface r\ =constant, 
and i — 1,2,3. In terms of 3-vectors, Eqns.Q-l^nj be- 
come 



V 2J D 4 = 0, 



V ? ;B 4 = 0, 



(11) 



(12) 



Advantage will be taken in the following of the formal 
equivalence of Eqns.J2J) and © and the equations for the 
electromagnetic field in a material medium, an equiva- 
lence that is fulfilled if we define 



= f{w)F- 



(4) 



A reference frame must be chosen in order to define the 
electric and magnetic fields from and F^ u . Since the 
cosmological model to be used in this article is described 
by the Friedmann-Robertson- Walker (FRW) metric, we 
can use the four velocity vector field = 6^ , which 
is orthogonal to the three-dimensional surfaces of homo- 
geneity of the FRW geometry. The electric and magnetic 
parts of the electromagnetic tensor are defined as: 



da 



(B l Y - —B i - rf ln ViE n = 0, 
da 



(13) 



(14) 



where the prime denotes derivative wrt conformal time. 

Since we are dealing with vector quantities we will use 
a vector basis P l defined by the following relations [2^| 

pi = pi( a; i) ) 



7 ZJ ViVjf* = -m l P\ 



D a = pa/3 y 



H a = P af3 V 01 



(5) 



7 y ViPj = 



E a = F afi V/3, B a = F af3 V/}. 



The eigenvalue m denotes the wave number of the corre- 
(6) sponding vector eigenfunction of the Laplacian operator 
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on the spatially homogeneous hypersurfaces. The spec- 
trum of eigenvalues depends on the 3-curvature e, and is 
given by 



rn 



2, < s < oo, e 



-1, 



^(1) = #S(1). (22) 
In terms of these, Eans. (|19|l and (|20|l can be written as 

P%' + (rn* + 2e)f(u)q%=Q, (23) 



m = s, < s < oo, e = 0, 

m 2 = s 2 - 2, s = 2,3, .., e = 1. 

The pseudovector basis J", convenient to work with the 
magnetic field, is defined by: 



(15) 



In terms of the basis P % and of the associated basis P l , 
the electric and magnetic fields can be expanded as fol- 
lows: 



l,m,(T 



(16) 



(17) 



These expressions are valid in the case e = 1. For e = 
and —1, the sum must be replaced by an integration 
((27r)~ 3 / d 3 x). The index I = 1,2 describes the two 
transverse degrees of freedom of the electric and mag- 
netic fields. Since we are using an expansion in standing 
waves |23, the index a takes the values "+" or "— ". 
Consequently, the fields in the case e = can be written 
as 

%) = Y1 ( E ™hil) cos(m . x) + E { J{ri) sin(m . x)) e|, 

(18) 

(where the e l l are modulus-one polarization vectors), and 
an analogous expression for B 1 . In the more general case 
of e = 1,-1 the (considerably more involved) analytic 
expression for the vector base can be found in [27] . The 
result of substituting in Eans. i|TT|) - ((rl|) the expansions 
given in Eans. l|16[l and l|17l) is 



*£ ) ' + T*& ) -/" 1 («)i>& ) =0 > 



(19) 



D%' + 7 D<$ + (m 2 + 2e) f(u) flg = 0. (20) 



The system described by Eans. l|19|l and i|20[l is not Hamil- 
tonian when written in the variables D and B, but a 
Hamiltonian can be introduced in terms of a new set of 
variables defined by 



(21) 



(a)' 
1ml 



(24) 



Since the new variables constitute a pair of canonical 
variables, the dynamical system given by Eans. p9(l and 
(I20|l can now be described by the Hamiltonian 



1 (<t)2 



61 9, 



(<7)2 

mi ' 



(25) 



where 



h(v)=f(co) 



6 2 (r?) = (m 2 + 2e)/(w). (26) 



The Hamiltonian for the mode characterized by (a, m, I) 
given in Eqn. (|25[l is identical to that of a single harmonic 
oscillator problem with time-dependent coefficients. We 
shall review in the next section how to quantize this sys- 
tem. 



II. QUANTIZATION 

The expression for TL obtained in the previous section 
appears in many different branches of physics, ranging 
from quantum optics to gravitational waves [28j. In or- 
der to quantize it, we shall follow a standard procedure 
of quantum optics. First we define creation and annihi- 
lation operators by the expressions 



Pj = ~ % 



(°0 



l u ml 



4<r) 
1ml 



+ < 



27 



(27) 



(28) 



where 7 is a constant to be determined later. In terms 
of d and d^ , the Hamiltonian operator can be written as 



1 



1 Hv) (2iV, 



where 
k(rj) = 761(17) 



+ l) + 7%) 4? + 7M1) 



;Wt2 

(29) 



62 (r?) 



62 fa) 



761(77), 



and M7 = d 



1/ 



is the number operator. The first 
term in Ean. l|29|l looks like the free Hamiltonian for the 
mode (cr, m, I) with a time-dependent function instead of 
the frequency. This part conserves the number of pho- 
tons in each mode, hence the total number of photons 
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and the total energy. The second and third terms rep- 
resent instead a time-dependent interaction which does 
not conserve the photon number. 

The operators aJ , and a m i obey the commutation re- 
lations 



a J. > a 



( cr )'tl _ r c 



for e = 1, 



(30) 



<5;z' S aa 'S(m — m ), fore = 0,-1, (31) 



where we have settled 7 = to in order to have vacuum as 
initial state. Since we shall work with a single mode (and 
in order to avoid clumsy expressions) the indices to, I and 
a will be omitted from now onwards. 

We will proceed to solve Schrodinger's equation for the 
Hamiltonian given in Eqn.J2HJ) by writing the time evolu- 
tion operator as a product of the rotation and the single- 
mode squeezed operators, along with a phase factor: 



U(r), 7/0) = e*S(ri)R{v), 
where S(rj) and R(i]) are defined as 



R{rf) 



ir(j7)aT 



(32) 



(33) 



(34) 



The function r(7y) is called the rotation angle, and ^(77) 
is usually defined as 



A{rf) = — r(rj) e 



(35) 



where r is the squeeze factor (defined in the range < 
r < 00) and tp the squeeze angle (defined in the range 
— 7r/2 < ip < 7r/2). The parameter r determines the 
strength of the squeezing while ip gives the distribution of 
the squeezing between conjugate variables. Schrodinger's 
equation with the Hamiltonian given in Ean. (|29|l reduces, 
after a direct but somewhat long calculation, to the fol- 
lowing system of first order coupled differential equations 
valid for each mode of the field: 



r' + (ip' + T') sinh2r = e 2iip 



ir' + (<p' + T') sinh2r = e 
(^'+r') cosh2r -ip' = 

V 2 



2' 



(36) 



(37) 



(38) 



(39) 



Similar equations were obtained by Albrecht et al 29] in 
the case of the scalar degrees of freedom of metric pertur- 
bations in cosmology. Grishchuk [3lT | derived analogous 



equations for the case of gravitational waves, and Mat- 
acz [3^, for a scalar field. To solve the system, we shall 
follow the method presented in . With the change of 
variables 



^(77) = e coshr(r7), 



(40) 
(41) 



P(rj) = _ e -%W+r/J) sinhr(77), 
Eans. H36|l - H39|l can be written as 

2a' = -ih/3 - ika, (42) 



2(3' = ikj3 + iha. 
Introducing the function fi(r]) defined by 

P* - a* 



13* 



Eans. H42|l and l|43|l are equivalent to 

2fi' - i{k - h)^ 2 + i(k + h) = 0. 
Finally, setting 



lh g 



(43) 



(44) 



(45) 



(46) 



we obtain the equation 



g" - %' + 6162.9 - 0. 

The problem has been reduced then to solving a single 
second-order differential equation for the function g{rf) 1 
which in terms of f(u>) reads (see Ean lj26|l ) 

g" + j^u J 'g' + (m 2 + 2e)g = 0. (47) 

We can see from this equation that there is no photon 
creation if the scalar field is constant, a fact that jus- 
tifies the identification of ui with the laser "pump" used 
in experiments devised to observe squeezed states [28^ 
With the function g we can get the squeeze parameter r 
through Eqn.gni as follows. From Eqns ijiOjl. JH}, and 
iHD we get 



l + ^W tanh r(r?) 



-l + e?Mv) tanh r (77)' 
which can be inverted to get r, given by 

i + M»?) + /**(»/) + IM»?)I 2 



tanh r(rj) 



1 



(48) 



(49) 



■ fi(T)) - fi*(rj) + \v{v)\ 2 ' 

All the functions appearing in the evolution operator can 
be written in terms of \i (we refer the reader to for 
details). A quantity of interest in the following is the 
mean number of photons per mode as a function of the 
conformal time, given by |35l | 



(N(r))) = sinh 2 r(?7). 



(50) 



In the next sections we shall calculate (N) for two differ- 
ent models. 



5 



III. A SIMPLIFIED MODEL 

Before going into the details of the model which takes 
into account the effect of matter creation on the metric, 
we shall study the production of photons in a simpler 
case, namely that in which the sole matter content is 
given by the scalar field u>. As discussed in there are 
nonsingular solutions in the theory of WIST that describe 
a FRW geometry plus a Weyl scalar field. Nonsingular 
solutions are also present in string theory (34|. Let us 
briefly review these solutions. The EOM for gravitation 
plus scalar field written in conformal time are |33| 



\ 2 

'2 i 2 , A I / \2 n 

a + ta + — [w a) = 0, 

D 



oj 1 = aa 2 , 



(51) 



where a = constant, and A 2 is the coupling constant of 
the scalar field to gravity. Notice that the equation of 
state of the scalar field is given by — , where 




FIG. 1: Plot of the scale factor (full line) and u as a function 
of the conformal time for A = |&| = 1. 



A 2 

From these equations we get 



9 ' 



(52) 



where we have defined b 2 — \ 2 a 2 /6. Eqn.JS^J) shows that 
only solutions with e = — 1 arc possible. Hence, from now 
on we shall restrict to the negative curvature case, for 
which Ean. (|52|l can be easily integrated. Although this 
choice is at odds with observation, we emphasize that 
we are using this exact solution only as a toy model to 
study the effect of matter creation on the production of 
photons. The result of the integration for the scale factor 
is 39], 



a(v) = V\bWcosh(2 V + 6). 



(53) 




From Ean. lfBT)! . 



FIG. 2: Plot of the mean number of photons as a function of 
the conformal time for m = 20 for the case without matter 
(dashed line) and for the case with matter creation (full line, 
see next section) , for A = 1 . 



ui(ri) — ±— — arctan (e 
2 A 



2r,+5\ , 7 J_ 



7T sqrt6 



(54) 



where 5 is an integration constant. The plots for these 
functions are given in Fig^ The scale factor displays 
a bounce, produced by the violation of the strong en- 
ergy condition by the scalar field |37^ . The number of 
photons created by the expansion in this model must 
be calculated numerically, since Ean. (|47|l has no analyt- 
ical solution for the scalar field given in Ean. (|54|l with 
the coupling f(w) = er 2uJ . The results are given by the 
dashed curve of Fig.J2l, for a wavenumber typical of the 
size of the intergalactic scale (1 Mpc). 



IV. A MODEL WITH MATTER CREATION 

We shall consider in this section the production of pho- 
tons in a background of the system composed of the scalar 
field plus matter and geometry, using a non-singular new 
solution that incorporates the effect of the creation of 
matter on the geometry. Friedmann's equation in con- 
formal time for this case is given by 



— [iv a) +—p m , 



(55) 
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while the second Einstein equation is 

( a" a 12 1 \ 

- 3 2— T y = p m + 3puj. (56) 

\ a J a* a -2 / 

The conservation of the stress-energy tensor gives 

d da 3 
-r(a 3 {p m +p u )) + (pm +Pw)~r- = o. (57) 
c(77 dry 

In the case of ultrarelativistic matter, this equation takes 
the form 

( fl W + ^(aV)'=0. (58) 

We would like to have a solution that describes creation 
of relativistic matter only around the bounce, and en- 
ters a radiation phase with a constant scalar field in a 
short time. Clearly, an asymmetry is to be expected both 
in the scale factor and in u>, since the evolution of this 
universe starts from the vacuum and enters a radiation- 
dominated epoch. An expression for a that fulfills these 
requirements is 

0(17) = /3v/cosh(2ry) + k sinh(2?y) - 2fc (tanh(r/) + 1), 

(59) 

with j3 = a v\J j^j^; an d < fco < 1/7- The expression 
for uj can be obtained using this scale factor in Ecins. l55|) 
and We shall not give the explicit expression, but 

the plots for 0(77) and w(r?) (see Fig|3J). Notice that the 
plot shows the announced asymmetry. The evolution of 
the Hubble parameter for the vacuum case (studied in 
the previous section) and for the case with matter cre- 
ation is showed in FigQ] Let us emphasize that since the 



scalar field tends rapidly to a constant value, the pro- 
duction of matter (controlled by uj', see Ecm.lj5%)l) stops 
soon, and the model enters a radiation phase without the 
need of a potential. In this sense, this solution describes a 
hot bounce, as opposed to cold bouncing solutions, which 
do not enter the radiation era unless they are heated up 
[3(| . Another nice feature of this solution is that the dila- 
ton goes automatically to a constant value for r\ — > 00, 
in such a way that the solution could be taken as the 
leading order of a perturbative development (as is the 
case in string theory). Again, no potential was needed 
in order to display this feature. The mean number of 
photons per mode can be calculated as in the previous 
cases, after numerical integration of Ean. (|47J) . The full 
line plot in Fig[3 shows the that the production of pho- 
tons is increased in the case of the model with matter 
creation. The number of photons is directly related to 
the fraction of electromagnetic energy stored per mode, 
relative to the background radiation energy /? 7 through 
the expression 

r(i/) « 

Pi 

where v = \Jm 2 — 2/a. Consequently, the increment in 
the number of photons may be of importance in the prob- 
lem of the generation of magnetic field seeds. 



V. DISCUSSION 

Using a covariant description for the electromagnetic 
field coupled to a scalar in a curved background, we ob- 
tained the Hamiltonian for each mode of the field. With 
the aid of the formalism of the squeezed states, we have 
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calculated the number of photons in a non-singular uni- 
verse using first an already known solution |33( with no 
matter present, and afterwards a new solution which 
takes into account the effect of the creation of ultrarel- 
ativistic matter on the evolution of the scale factor and 
of the scalar field. This new solution presents several in- 
teresting features, namely the transition of an expanding 
phase to a radiation era, and the constancy of the scalar 
field a short time after the bounce. The results for the 
mean photon number per mode show that the produc- 
tion is increased in the case of matter production. This 



increment may be relevant for the creation of seeds of the 
magnetic field. We hope to discuss in detail this issue in 
a forthcoming publication. 



Acknowledgements 

SEPB acknowledges UERJ for financial support. JS is 
supported by CNPq. 



[1] E. Parker, Ap. J. 162, 665 (1970). See also Y. B. Zel- 
dovich, A. A. Ruzmaikin, and D. D. Sokoloff, Magnetic 
Fields in Astrophysics, Mc Graw-Hill, New York (1980). 

[2] See however R. Kulsrud, S. C. Cowley, A. V. Gruzinov 
and R. N. Sudan, Phys. Rep. 283, 213(1997), and L. M. 
Widrow, Rev. Mod. Phys. 74, 775 (2003). 

[3] L. Biermann, Z. Naturf. 5A, 65 (1950). 

[4] T.E. Clarke, P.P. Kronberg and H. Bohringer, Astrophys. 
J. 547, Llll (2001). 

[5] H. Bohringer, Rev. Mod. Astron. 8, 295 (1995). 

[6] S. A. Colgate and H. Li, |astro-ph/0001418| 

[7] D. Grasso and H. Rubinstein, Phys. Rept. 348, 163 
(2001). 

[8] M. Turner and L. Widrow, Phys. Rev. D 37, 2743 (1988). 
[9] For a review on different mechanisms of seed generation, 

see Generation of magnetic fields in cosmology, A. Dol- 

gov, in Gurzadyan, V.G. (ed.) et al.: From integrable 

models to gauge theories, 143-154, hep-ph/01 10293. 
[10] Gauge invariant wave equations in curved space-times 

and primordial magnetic fields, G. Lambiase and A. 

Prasanna, gr-qc/0407071. 
[11] M. Novello and J. Salim, Phys. Rev. D 20, 377 (1978). 
[12] A. Dolgov, Phys. Rev. D 48, 2499 (1993). 
[13] E. Calzetta, A. Kandus and F. Mazzitelli, Phys. Rev. D 

57, 7139 (1998). 
[14] B. Ratra, Astrophys. J. 391, LI (1992). 
[15] Nearly minimal magnetogenesis, T. Prokopec and Ewald 

Puchwein, astro-ph/ 0403335. 
[16] Lemoine D., Lemoine M., Phys. Rev. D 52, (1995), 1955, 

Gasperini M., Giovannini M. and Veneziano G. Phys. 

Rev. Lett. V 75, (1995), 3796. 
[17] J. Lidsey, D. Wands, E. Copeland, Phys. Rept. 337, 343 

(2000), hep-th/9909061. 
[18] J. M. Salim and S. Sautu, Class. Quant. Gravity, 13 

(1996) 353-360, J. M. Salim, S. Sautu and Martins R., 

Class. Quant. Gravity, 15 (1998) 1521. 
[19] P. Dirac, Nature 139, 323 (1937). 

[20] See Magnetogenesis, variation of gauge couplings, and 
inflation, M. Giovannini, Proceedings of the Chalongc 
School on Astrofundamental Physics, N.G. Sanchez and 



Y.M. Pariiski, eds. Kluwer (2002), astro-ph/0212346. 

[21] K. Bamba and J. Yokoyama, Phys. Rev. D 69, 043507 
(2004), astro-ph/0310824. 

[22] Large-scale magnetic fields from dilaton inflation in non- 
commutative spacetime, K. Bamba and J. Yokoyama, 
hep-ph/0409237. See also Generation of Cosmologi- 
cal seed magnetic fields from inflation with cutoff, A. 
Ashoorioon and R. Mann, gr-qc/0410053. 

[23] See for instance D. Bailin and A. Love, Rept. Prog. Phys. 
50, 1087 (1987). 

[24] See for instance T. Damour and A. M. Polyakov, Nucl. 
Phys. B423, 532 (1994), hep-th/9401069. 

[25] Lifshitz E. M. and Khalatnikov I. M., Adv. Phys. 12, 
(1963), 185. 

[26] M. Novello, J. M. Salim, M.C. Motta da Silva, S. E. Joras 
and R. Klippert, Phys. Rev. D 52, (1995), 730. 

[27] E. Kalnins and W. Miller, Jr., J. Math. Phys 32, 698 
(1991). 

[28] B. Schumaker, Phys. Rep. 135, 317 (1986). 

[29] A. Albrecht, P. Ferreira, M. Joyce, T. Prokopec, Phys. 

Rev. D 50, 4807 (1994), astro-ph/9303001 
[30] L. Grischuk and Y. Sidorov, Phys. Rev. D 42, 3413 

(1990). 

[31] L. Grishchuk, Phys. Rev. D 45, 4717 (1992). 

[32] Matacz A. L., Phys. Rev. D 49, (1994), 788. 

[33] M. Novello, L. Oliveira, J.M. Salim and E. ELbaz, Int. 
J. of Mod. Phys. D 1, 641 (1993). 

[34] See for instance Old ideas and new twists in string cos- 
mology, M. Giovaninni, hep-th/0409251. 

[35] M. Giovannini, Primordial Magnetic Fields, 
hep-ph/0208152. 

[36] M. Giovannini, Class. Quantum Grav.21, 4209 (2004). 

[37] C. Molina-Paris and M. Visser Phys. Lett. B455, 90 
(1999). 

[38] Notice that the expansion in standing waves is completely 
equivalent to that of travelling waves (see for instance 

E3). 

[39] The expression for the scale factor in terms of cosmolog- 
ical time was given in js^ ]. although not in a closed form 



